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Recently a holographic model of hadrons motivated by AdS/CFT has been proposed 
to fit the low energy data of mesons. We point out that the infrared physics can be 
developed in a more systematic manner by exploiting backreaction of the nonperturba- 
tive condensates. We show that these condensates can naturally provide the IR cutoff 



Oh 

| corresponding to confinement, thus removing some of the ambiguities from the original 

43 ■ formulation of the model. We also show how asymptotic freedom can be incorporated 

into the theory, and the substantial effect it has on the glueball spectrum and gluon 
condensate of the theory. A simple reinterpretation of the holographic scale results in a 
non-perturbative running for a s which remains finite for all energies. We also find the 



leading effects of adding the higher condensate into the theory. The difficulties for such 
models to reproduce the proper Regge physics lead us to speculate about extensions 
of our model incorporating tachyon condensation. 



1 Introduction 



QCD is a perennially problematic theory. Despite its decades of experimental support, the 
detailed low-energy physics remains beyond our calculational reach. The lattice provides a 
technique for answering nonperturbative questions, but to date there remain open questions 
that have not been answered. For instance, the low-energy scalar spectrum is a puzzle. 
There are a lot of experimentally observed states, however their composition (glueball vs. 
quarkonium) and their mixings are not well understood. The difficulty for any theory try- 
ing to make progress in this direction is to understand the interaction between the scalar 
states and the vacuum condensates of QCD. In this paper we attempt to incorporate the 
effects of the vacuum condensates into the holographic model of QCD as a first step toward 
understanding the scalar sector in the context of these models. 

The SVZ sum rules [1] are a powerful theoretical tool for relating theoretically solid 
facts about perturbative QCD with experimental data in the low-energy region. The basic 
observation is that the correlator for a current J 



may be expanded in a Wilson OPE that is valid up to some power (where instanton cor- 
rections begin to invalidate the local expansion), (l/Q 2 )(dn/d(5 2 ) = ^C2d(02d), where 
the Oid are gauge-invariant operators of dimension 2d. In the deep Euclidean domain the 
coefficients are calculable. On the other hand, the correlator relates to observable quan- 
tities; for instance if J is the vector current then Imll(s) is directly proportional to the 
spectral density p(s), measurable from the total cross section a(e + e~ — > hadrons). The SVZ 
sum rules, then, relate Wilson OPE expansions to measurable quantities. They can become 
useful for understanding detailed properties of the first resonances when one takes a Borel 
transform that suppresses the effect of higher resonances. As it turns out, keeping only low 
orders of perturbation theory in the coefficients C^d, one still obtains reasonable agreement 
with data, so the assumption that the largest corrections arise from the condensates is a 
good one. 

Recently another technique for understanding the properties of low-lying mesons has 
arisen in the form of AdS/QCD. The phenomenological model constructed on this basis [2] 
takes as its starting points the OPE (much as in the SVZ sum rules) and the AdS/CFT 
correspondence [3]. The idea is straightforward: rather than attempting to deform the 
usual Type IIB on AdS$ x S 5 to obtain a theory more like QCD, one starts with QCD 
and attempts to build a holographic dual. Of course in detail such a program is bound 
to eventually run up against difficulties from a' corrections, g s corrections, the geometry 
of the five compact dimensions (or the proper definition of a noncritical string theory) and 
other issues. However, one can set aside these problems, begin with a relatively small set 
of fields needed to model the low-lying states in QCD, and see how well the approximation 
works. In this phenomenological approach with bulk fields placed in the Randall-Sundrum 
background [4], and the AdS space cut off with a brane at a fixed z = z c in the infrared 
a surprisingly good agreement with the physics of the pions, p, and a± mesons has been 
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found [2] (for more recent work on AdS/QCD see [5,6]). However, there are several obvious 
limitations to this model. For example, it does not take into account the power corrections 
to the OPE in the UV (the effects of the vacuum condensates), or the corrections coming 
from logarithmic running. Also, the theory has a single mass scale (set by the location of IR 
cutoff brane) which determines the mass scales in all the different sectors in QCD. 

Here we show how the effect of the vacuum condensates and of asymptotic freedom can 
be simply incorporated into the model. We will limit ourselves to pure Yang-Mills theory 
without quarks, though we expect that most of the aspects of incorporating quarks should 
be relatively straightforward. For the vacuum condensates one has to introduce a dynamical 
scalar with appropriate mass term and potential coupled to gravity. There will be a separate 
field for every gauge invariant operator of QCD, and a non-zero condensate will lead to a 
non-trivial profile of the scalar in the bulk. While the effects of these condensates on the 
background close to the UV boundary are small (though not always negligible), they will 
become the dominant source in the IR and effectively shut off the space in a singularity (and 
without having to cut the space off by hand). This resolves several ambiguities in choosing 
the boundary conditions at the IR brane, and the various condensates will also be able to set 
different mass scales in the different sectors of QCD. Asymptotic freedom can be achieved 
by properly choosing the potential for the scalar corresponding to the Yang-Mills gauge 
coupling. Incorporating asymptotic freedom will have an important effect on the glueball 
sector since a massless Goldstone field will pick up a mass from the anomalous breaking of 
scale invariance. 

It is natural to ask why we should expect QCD to have any useful holographic dual. The 
most well understood examples of holography are in the limit of large N and gy M iV ^> 1, 
far from the regime of real-world QCD, which apparently would be a completely intractable 
string theory with a large value of a'. On the other hand, the large N approximation has 
frequently been applied in QCD phenomenology. A major concern in applying holography 
to QCD is whether the dual should be local, or whether it can have higher-derivative a' 
corrections. The a' corrections are associated with the massive stringy excitations in the 
bulk, which once integrated out yield complicated Lagrangians for the remaining light fields. 
However, we can think of this physics in a different way. Holography is closely related to 
the renormalization group [7]; the coordinate z can be identified with where /i is the 
renormalization group scale. AdS/CFT identifies massive bulk fields with higher-dimension 
operators in the field theory. From this point of view, a' corrections to the physics of light 
bulk fields are associated with the effects of higher- dimension operators coupled to the low- 
dimension operators in the RG flow. The OPE tells us that, at large Q 2 , the effects of these 
higher-dimension operators in the field theory are controllably small. From this point of 
view, despite the apparently large a' corrections, it is reasonable to begin with a local bulk 
action in terms of fields corresponding to the low dimension operators of QCD. The example 
of SVZ gives us hope that this can correctly capture physics of light mesons. For highly 
excited mesons, which in QCD look like extended flux tubes, and thus feel long distances, it 
is more probable that the strong IR physics will mix different operators and that our neglect 
of a'-like corrections will become more troublesome. It is a general problem, in fact, that 
physics of highly excited hadrons is troublesome in these models [8]. Recent proposals for 
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backgrounds with correct Regge physics [9] offer some hope of addressing this problem. We 
will offer some further comments on how a closed-string tachyon might give a dynamical 
explanation of such a background, but there is clearly much more to be done along these 
lines. 

The paper is organized as follows. In section 2 we will remind the reader of the basic 
formulation of AdS/QCD on Randall-Sundrum backgrounds as in Ref. [2], and discuss some 
of its shortcomings. In section 3 we will discuss a class of models that do not incorporate the 
running of the QCD coupling, but do incorporate the effects of the lowest vacuum condensate. 
We will point out that these backgrounds too have some shortcomings, but improve on the 
RS backgrounds, and may provide a useful setting for exploring some questions. We calculate 
the gluon condensate and point out that there is a zero mode in the glueball spectrum due 
to the spontaneous breaking of scale invariance. In section 4 we discuss the construction of 
5D theories with asymptotic freedom. We calculate the gluon condensate and provide a first 
estimate of masses for the ++ glueballs from these backgrounds. We point out a possible 
interpretation of the running coupling where a s remains finite for all energies. In section 5 
we show how to systematically incorporate the effects of the higher condensates. We give a 
background including the effect of (TrF 3 ) and show how it affects the gluon condensate and 
the glueball spectrum. In section 6 we discuss the difficulties of reproducing the correct Regge 
physics, and speculate that closed string tachyon dynamics could perhaps be responsible for 
reproducing the necessary IR background. Finally we provide the conclusions and some 
outlook in section 7. 



2 AdS/QCD on Randall-Sundrum backgrounds 

We briefly review the AdS/QCD model on Randall-Sundrum backgrounds [2]. One assumes 
that the metric is exactly AdS in a finite region z = to z = z c , i.e. 




where z c ~ A-qcd determines the scale of the mass spectrum. It is assumed that there is a 
brane (the "infrared brane") at z — z c ; in practice one assumes a UV boundary at z = e 
and sends e — » at the end of calculations. One puts a field <po in the bulk for every gauge 
invariant operator O in the gauge theory. If O is a p-form of dimension A, <fi has a 5D mass 
m\ = (A + p)(A + p — 4). This operator has UV boundary conditions dictated by the usual 
AdS/CFT correspondence. On the other hand, IR boundary conditions are less clear, and 
one in principle can add localized terms on the IR brane. 

For instance, in the original papers the treatment involves the rho and a,\ mesons and the 
pions. There are bulk gauge fields, Aj^ and (where M is a 5D Lorentz index) coupling 
to the operators <?L7^ a <?L and q~R^t a qR. There is also a bulk scalar X alS coupling to the 
operator q^q^, where a and (3 are flavor indices. The scalar X alS is assumed to have a 
profile proportional to 5 al3 (m a z + (q a q a ) z 3 ), based on the Klebanov-Witten result [11] that 
a nonnormalizable term in a scalar profile corresponds to a perturbation of the Lagrangian 
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by a relevant operator (in this case, mqq), while a normalizable term corresponds to a 
spontaneously generated VEV for the corresponding field theory operator. The profile for 
X(z) couples differently to the vector and axial vector mesons and achieves the p — a\ 
mass splitting. Furthermore, the pions arise as pseudo-Nambu-Goldstone modes of broken 
chiral symmetry, and the Gell-Mann-Oakes-Renner relationship is satisfied. Several other 
constants in the chiral Lagrangian are determined to around 10% accuracy. 

There are several limitations to this model. One is that it does not take into account 
the power corrections to the OPE in the UV, or the corrections coming from logarithmic 
running. The SVZ sum rules work fairly well with just the leading power correction taken 
into account [1], and these power corrections can be incorporated into the form of the metric 
with a simple ansatz [6]. However, backreaction has not been taken into account in such 
studies, so the Einstein equations will not be exactly satisfied. Also, in the SVZ sum rule 
approach, leading corrections to the OPE essentially determine the mass scale of the lightest 
resonance in the corresponding channel. In the Randall-Sundrum approach, the leading 
correction to the OPE and the IR wall at z c both influence the corresponding mass scale. 
This has effects on the spectrum. 

For instance, in QCD the mass scales associated with mesons made from quarks and 
mesons made from gluons are very different [10]. One can see this in the OPE. For instance, 
in the tensor 2 ++ channel, the leading corrections to the OPE come from the same operator 
in the qq case as in the G 2 case, but the coefficients are very different. In the Randall- 
Sundrum approach, without turning on a background field VEV the mass scale for both of 
these channels is set by z~ x . One can turn on a background VEV and give the fields different 
couplings to it to reproduce the difference in the OPEs, but z~ 1 will continue to play a role 
in setting their masses. 

Another limitation is that these backgrounds lack asymptotic freedom. If one defines 
the theory with a cutoff at which a s has some finite value, this might not appear to be of 
central importance. On the other hand, we understand that in real QCD the values of the 
condensates are determined by the QCD scale at which the perturbative running coupling 
blows up, so achieving asymptotic freedom can allow such a relationship and make the model 
less ad hoc. Furthermore, we will see that the lack of proper conformal symmetry breaking 
can lead to a massless scalar glueball state (i.e., the model has a "radion problem"), which 
is most satisfactorily resolved by incorporating asymptotic freedom. 

In summary, the AdS/QCD models on hard- wall backgrounds work surprisingly well 
for some quantities, but have obvious drawbacks. There are ambiguities in IR boundary 
conditions, and the existence of a single IR wall influencing all fields obscures the relationship 
between masses of light resonances and power corrections discovered by Shifman, Vainshtein, 
and Zakharov. Luckily, there is a simple remedy to these difficulties: we simply remove the 
IR brane, and allow the growth of the condensates to dynamically cut the space off in the 
infrared. 
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3 Vacuum condensates as IR cutoff 



To model the pure gauge theory we begin with the action for five-dimensional gravity coupled 
to a dilaton (in the Einstein frame): 

1 C 19 1 

S = ^ J d 5 x^(-U + + l -g MN d M <t>d N <t>). (3.1) 

Here k 2 is the 5 dimensional Newton constant and R is the AdS curvature (related to the the 
(negative) bulk cosmological constant as R~ 2 = — ^-A. Note, that <fi is dimensionless here. 
The dilaton will couple to the gluon operator G lw G lw '. The fact that there is a non- vanishing 
gluon condensate in QCD is expressed by the fact that the dilaton will have a non-trivial 
background. We can find the most general such background by solving the coupled system 
of the dilaton equation of motion and the Einstein's equation, under the ansatz that we 
preserve four-dimensional Lorentz invariance while the fifth dimension has a warp factor: 

ds 2 = e-^rj^dx^dx" - dy 2 

= 0(2/)- (3.2) 
The coupled equations for A(y), <f>(y) will then be 

AA' 2 - A" = AR 2 

A ~ 24 + R 2 

0" = AA'(j)'. (3.3) 

A simple way of solving these equations is to use the superpotential method [12, 13], that is 
define the function W{4>) such that 

A'(y) = W{4>{y)) 
dW 

#{V) = 6^- (3-4) 
This is always possible, and the equation determining the superpotential is given by 

dW\ 2 o 12 



To solve for the most general superpotential consistent with our chosen (constant) potential, 
it is useful to parameterize it in terms of a "prepotential" w as [13] 

w = 4 f w+ 1 



R\ w 

W>=ML-l), (3.6) 
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which is chosen such that ()3.5|) is automatically satisfied. The consistency condition of the 
two equations in ()3.6|) implies a simple equation for the prepotential: 



w 



w. 



and thus for the superpotential we find 



R 



(3.7) 



(3.8) 



With the superpotential uniquely determined (up to a constant C) we can then go ahead 
and integrate the equations in (|3.4|) . The result is 



<Ky) = y-log[Ctanh2(y -y)/i?], 

A{y) = log [cosh2(y -y)/Rsmh2{yo-y)/R]+A , 



(3.9) 



where y$ and A are integration constants. These are the solutions also found in [13,14]. 
In order to have A(y) asymptotically equal to y (for large negative y), we will fix A = 
^ — ~ log 2. Note, that there is another branch of the solution where in which tanh is 
replaced by coth in the solution for 0. In order to have a form of the solution that is more 



familiar and useful, we make a coordinate transformation e h = z/z c . This will recast the 
solution in a form that matches the usual conformal coordinates z near z = 0: 



ds 2 




1 - ( — ) r ]ilv dx il dx v - dz< 



1 - 



+ 



(3.10) 



(3.11) 



/ 



Here z c is a new parameter determining the IR scale, and we expect z c ~ ^qcd- The 
point z = z c is a naked singularity, which we must imagine is resolved in the full string 
theory. Also note that the first correction to the AdS 5 metric goes as z 8 . 



3.1 Gluon condensate 

We have seen that the metric incorporates power corrections to the pure AdS solution, which 
we want to identify with the effects of the gluon condensate. Below we would like to make 
this statement more precise. According to the general rules of the AdS / CFT correspondence 
given some field (po on this background representing an operator O in QCD, the deep- 
Euclidean correlator of O will have a Q~ 8 correction if ipo has no dilaton coupling, and a 
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Q 4 correction if ipo couples to 0. Note that near z = 0, the dilaton behaves + This 
is in agreement with the expectation that a field coupling to an operator with dimension A 
has two solutions, z A ~ d and z A . In our case A = d — 4, and we expect the constant piece 
to correspond to the source for the operator TrG 2 and coefficient of the z 4 to give the gluon 
condensate. The precise statement [11] from AdS/CFT is that if solution to the classical 
equations of motion $ has the form near the boundary 

z) -> z d - A [%(x) + 0(z 2 )] + z A [A(x) + 0(z 2 )] (3.12) 

then the condensate (one-point function) of the operator O is given by 

(0{x)) = (2A -d)A(x). (3.13) 

However, to apply this to the solution in (|3.11j) we need to make sure that the appropriate 
normalization of the fields is used. The expression ()3.13)1 is derived assuming a scalar field 
action of the form 1/2 J d A xdz / V ! (9jy</>) 5 ! . Comparing this with the action used here (|3.1|) 
we find 

(TrG 2 ) = V3a/ — -. (3.14) 

V K Zl 

In order to be able to relate this expression for the condensate we need to find an expres- 
sion for R 3 / k 2 . This can be done by requiring that the leading term in the OPE of the gluon 
operator G 2 is correctly reproduced in the holographic theory. One can simply calculate the 
leading term in the OPE in QCD [15] 

J d 4 x(G 2 (x)G 2 (0)}e^ = - (iV ^~ 1) g 4 log^ + . . . (3.15) 

The same quantity can be calculated in the gravity theory by evaluating the action for the 
scalar field with a give source <fio(q) in the UV. The general expression for the action is 
obtained (after integrating by parts and using the bulk equation of motion) 

If R 3 1 

S 5 d = 7T^ / d x^r-(j)d z (j)\ z ^ Q . (3.16) 

ZKT J Z 6 2 

In order to find the action one can use the bulk equation of motion close to the UV for the 
scalar field given by 

<p" + q 2 (p - V = 0. (3.17) 

Requiring that the wave function approaches 1 around z = will fix the leading terms in 
the wave function: 

<f>{z) = 1 - igV log q 2 z 2 + l -q 2 z 2 + .... (3.18) 
Taking the second derivative of the action we find 

( d 4 x(G 2 (x)G 2 (0))e iqx = --^g 4 logg 2 //i 2 + . . . , (3.19) 
J lofir 
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from which we get the identification 



R3 - ^ ~ 1 >, (3.20) 



K 2 7T 2 



Using this result we get a prediction for the gluon condensate 

(TrG 2 ) = -^v/3(iV 2 - 1). (3.21) 



3.2 The Glueball Spectrum 

Now we wish to solve for the scalar glueball spectrum, which is associated with scalar fluctu- 
ations about the dilaton-metric background. The analogous calculation in the supergravity 
model has been performed in [16]. We should solve the coupled radion-dilaton equations, as 
often there is a light mode from the radion. In other words, we should solve for eigenmodes 
of the coupled Einstein-scalar system. This has been worked out in detail for a generic scalar 
background in [17]: the linearized metric and scalar ansatz is given by: 

ds 2 = e- 2A{y) (l - 2F(x,y))dx^dx u r ]llu - (1 + AF(x, y))dy 2 

<j>(x, y) = Mv) + 4- {F'(x, y) - 2A'(y)F(x, y)) . (3.22) 

K 2 (f) 

This will satisfy the coupled Einstein-scalar equations if F — F(y)e iq x with q 2 = m 2 and 
F(y) satisfies the differential equation 



F" - 2A'F' - 4A"F - 2^-F' + 4A'^F = -e" 2A m 2 F. (3.23) 

Using the solutions for A(y) and 4>o(y) from Eq. 13.91 and an ansatz F(x,y) = F(y)e iq ' x , with 
q 2 = —m 2 , this becomes: 

F"{y) + ^cothfcMF'(y) + f § + m 2 ^£^== | F(y) = 0. (3.24) 

R R V v 2sinhM ^. 

Demanding a normalizable solution, we need that (in the z coordinates) J dz^fg\(p(z)\ 2 and 
J dzy/gg 55 \d z (p(z)\ 2 be finite. Thus we need (p(z) ~ z 4 at small z. We need to be somewhat 
more careful about solving for F: in the z coordinates, we have that 

dF dA „ dF „ . , . . 

z^0:z- 2z—F — > z— 2F ~ <p(z), (3.25) 

dz dz dz 

so that F(z) ~ z 2 near z = is compatible with our assumptions on the behavior of (p(z). 
As usual, this equation can be solved using the shooting method: the differential equation 
is solved numerically starting from the UV boundary with arbitrary normalization the BC 
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following from (|3.25|) (the numerics is very insensitive to the choice of the actual UV BC for 
the higher modes) for varying values of m 2 . For discrete values of m 2 the wave function will 
be normalizable (which numerically is equivalent to requiring a Neumann BC at the location 
of the singularity). This way we find (in units of z~ l ) glueballs with masses 6.61, 9.84, 12.94, 
and 15.98 (and so forth, with regular spacing in mass). 

There is also a serious problem: we find a massless mode for the radion (with the F(z) ~ 
z 2 UV boundary conditions). This can be understood as follows: in real QCD, the classical 
conformal symmetry is broken by the scale anomaly. However, in our model, it is broken 
by the z 4 profile of the dilaton. AdS/CFT tells us we should understand turning on such 
a normalizable background in the UV as a spontaneous symmetry breaking. In Randall- 
Sundrum models, one similarly has a radion problem and needs to invoke (for instance) 
a Goldberger-Wise stabilization [18] to avoid a massless mode. The radion has not been 
part of previous investigations of glueballs on Randall-Sundrum backgrounds [19], so such 
studies have essentially assumed that the stabilization mechanism removes a light mode 
from the spectrum. However, an added Goldberger-Wise field does not seem to correspond 
to an operator of QCD, so it goes against the spirit of the AdS/CFT correspondence. It 
is apparent that a palatable solution of the radion problem in AdS/QCD demands a 5D 
treatment of scale dependence that mirrors that in 4D. This motivates us to search for 
backgrounds incorporating asymptotic freedom, which also allows us to begin approaching 
a more detailed matching to perturbative QCD. 

4 Incorporating Asymptotic freedom 

We have shown in the previous section a background that incorporates the lowest QCD 
condensate TrF 2 and which automatically provides an IR cutoff via the backreaction of the 
metric. However, this setup is certainly too simplistic even to just produce the main features 
of QCD: for example asymptotic freedom is not reproduced in that setup. It is the dilaton 
field that also sets the QCD coupling constant, and by approaching a constant value the 
model in the previous section actually describes a theory that approaches a conformal fixed 
point in the UV, rather than QCD. One may think that this is not an important difference 
for the IR physics, but this is not quite right. For example as we have seen it introduces a 
"radion problem." 

Thus we set out to find a potential for the dilaton that will reproduce the logarithmic 
running of the coupling. We assume that similarly to string theory the gauge coupling is 
actually given by e b ^ z ' (where b is a numerical constant). We will find a result consistent 
with expectations from string theory. 1 

We assume that our action is: 

S = ^ J d 5 x^(-K - V{4>) + ~£ MJV cM«), (4.1) 

where now we will try to determine V{<p) such that we reproduce asymptotic freedom. If we 
1 Other discussions of backgrounds with logarithmically running coupling can be found in [20,21]. 
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require that the coupling runs logarithmically, and as usual identify the energy scale with 
the inverse of the AdS coordinate z we need to have a solution of the form 

P b <t>{z) _ (A 9) 

6 "log 2»' 1 ] 

where z = Aqq D and we do not fix b a priori. Then going to the y coordinates as usual via 
the definition e y l R = z/R, we will find 

e b<f>(y) = _J_ ( 4 _ 3 ) 

yo-y 

If we now assume that this solution follows from a superpotential W then (j)'{y) = 6^ = 
b (yl_ y } — b^ e ^- This implies that W(<f>) = ^e^ + Wo. Now that we have found the form of 
the superpotential, we can easily solve for the warp factor: A'(y) = W(<p{y)) = 6b 2^_ yo ^ 
and hence A(y) = A + Way + gpTog ^z^- In z coordinates, this becomes: 

A(z) = A + W R\ogz/R--^\og\ogz /z, (4.4) 

and hence e~ 2A ^ = e~ 2A °(R/z) 2WoR (\ogZo/z) 1 ^ 3b2 \ From this we conclude that we should 
take A = and W = to get a solution that looks AdS-like up to some powers of log z Q /z. 
The potential corresponding to this superpotential is then given by 

™ = 18 (f ) 2 - mv ' 2 = -»h{(h-D e2W + 12eW + 36ft2 ) • (45) 

This is particularly simple in the case that b = ±-t/§- In that case we have simply 

V{4>) = ~ f (4-6) 



(4.7) 



= T^log^ = ^loglog^ (4.8) 

. y 1 . -R . z 1 , . , , _ N 

A=— + -log = log- - -log log— , 4.9 

R A y -y R 4 z 

and thus the metric will be 

rf s 2 = Q (^Jlogjdx^dx^ - dz 2 ^j = - j=^ dx»dx v r]^ - dy 2 . (4.10) 

Our dilaton in ()4.1j) is normalized in an unusual way, nevertheless (J4.6J) is recognizable 
as a potential that commonly occurs in nonsupersymmetric string theory backgrounds: a 
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cosmological constant plus a term exponential in the dilaton. This is quite reasonable from 
the string theory perspective, where such a term can arise from dilaton tadpoles in critical 
backgrounds or from the central charge in noncritical backgrounds. In fact, the factor a/2/3 
arises from string theory considerations in a simple way, which may be an amusing coinci- 
dence or may have more significance. Suppose that there is a noncritical string theory in 5 
dimensions. Its action in string frame has the form [22] 

S = -L f d 5 x(-G) 1/2 e- 2 * (C + R + Ad^d^ + •••), (4.H) 
2Wq J 

where C is proportional to the central charge and is nonvanishing since we are dealing with 
a noncritical string. Now we go to Einstein frame: 

S = J d 5 x(-G) 1/2 (ce^l* + R- p^d^ + ••-), (4.12) 
and finally we note that comparing to our normalization above, <£> = a/3/80, so that e 4 */ 3 = 



4.1 The Glueball Spectrum 

To calculate the glueball spectrum we can apply Eq. ()3.23|) for the background in tj4.8H4.9jl . 
This equation (transformed to z coordinates) reduces to (in units of R) 

^-m - * o + dy f '« + (i4 + vS?) fw = ° (4 - i3) 

Using the shooting method again we find (in units of z^ 1 ) glueballs at 2.52, 5.45, 8.16, and 
10.81. In particular this background seems to have a light mode from the radion, but not 
any zero mode. 

Lattice estimates put the first ++ glueball in pure SU(3) gauge theory at approximately 
1730 MeV, and the second at about 2670 MeV, with uncertainties of order 100 MeV [23]. 
Thus they put the ratio of the first and second scalar glueball masses about about 1.54, 
whereas we find a significantly larger value of 2.16. While the lattice errors are still fairly 
large, this probably indicates that we are not so successful at precisely determining properties 
of the second scalar glueball resonance. Since we undoubtedly fail to properly describe highly 
excited resonances, this is not so surprising. If we set z^ 1 to match the lattice estimate for 
the first glueball mass, we find 

z l » 680 MeV. (4.14) 

One can also calculate the spectrum of spin 2 ++ glueball masses by solving the fluctu- 
ations of the Einstein equation around the background. The resulting differential equation 
we find is 

zlog^f(z) - (l + 3\og^)f(z)+m 2 z x fk^f(z) = 0. (4.15) 
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Using the shooting method (and imposing Dirichlet BC on the UV boundary) we find the 
lightest modes at 4.03, 6.56, ... in units of l/z . It is important to point out that the lowest 
spin 2 ++ glueball is naturally heavier in this setup than the spin ++ glueball due to the 
mixing of the radion with the dilaton. In the usual supergravity solutions the spin ++ and 
2 ++ glueballs usually end up degenenerate (in contradiction to lattice simulations). This is 
for example the case in the AdS black hole solution of Witten analyzed in [16] (the additional 
light scalar modes identified in [24] do not correspond to QCD modes). 



4.2 Power Corrections and gluon condensate 

We would now like to evaluate the gluon condensate in this theory assuming that the the IR 
scale z is fixed by the value of the lightest glueball mass. As usual one needs to calculate 
the 5D action corresponding to a fixed source term turned on for the QCD coupling and to 
get the condensate (one -point function) we need to differentiate the 5D action with respect 
to the source. Ordinarily, the computation of the 5D action on a given solution in AdS/QCD 
reduces to simply a boundary term. However, in our case it is not so simple: our potential 
is not just a mass term, so that the bulk piece V — |^ is not set to zero by the equations 
of motion. As a result, the action also has a "bulk" piece. 

We evaluate the 5D action as a function of z , imposing a UV cutoff at e. The action is 
given by: 

i £ {li^l ("* - 5*™* + | + *. (4.16) 



This integral can be performed explicitly: 



1 

2^ 



1 2 log ^ 20 
I 

2z 4 z 



4 



(4.17) 



We drop the UV divergent terms (the 1/e 4 pieces) assuming that there will be counter terms 
absorbing these. Then the explicit expression for the action will be 

= (4.18) 

Note that this is more easily calculated as 

S(z ) = J^J d4x 2 V9lB W(<f>), (4.19) 

evaluated at the boundary z = z , where g±D is the induced 4D metric at the boundary. 
(This observation has been made before in Ref. [25].) To see this, we use the following 
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relations: 



dW 



2 



Tl(y) = -20 A' (y) 2 + 8 A" (y) = -20^ + 48 1 — 1 (4-20) 

->> 2 - -Kf ) 2 (4 ' 2i) 



-V{4>) = -18 ( — ) + 12W 2 (4.22) 
to see that i/gS^o, evaluated on the solution, is 



e 




-*A(v) ( _ 8 ^ 2 + 12 ( ^ ) ) = 2-| (e- a (%(^(|,))) . (4.23) 



This makes it clear that we can use the superpotential as a counterterm on the UV boundary 
to cancel the terms diverging as e — > (which we dropped above.) 

It turns out that \ is almost precisely as in the background without asymptotic freedom 
(because, for the fluctuating modes, corrections to the wavefunction near z = are small a s 
corrections), so one can still use ()3.20|) to find the value of R 3 /k 2 . However, there is a slight 
subtlety: we found a value for -=j assuming a source coupled to TrG 2 . In fact in our case 

we have fixed the numerical factor in the correspondence of eV^ 2 /^ based on its asymptotic 
behavior. Using the expression for the coupling in a pure YM theory 

3 Jv c iug Aqcd 

and identifying Aqcd = ^ and Q = -, we have at the cutoff z = e = 



f2,( A HN C , , , lliV r n . . 

= V ayM(A) = 24^^ (A) - (4 - 25) 



Now, for a fluctuation <p(z), we have 



2/3WW +!P W e ««-«) w eV /2/3^) (1 + y/^/Sytye**). (4.26) 

Now, <^(z) near z = behaves like any massless scalar fluctuation on an AdS background, 
and thus we have that it shifts the action by an amount 

Ssn = ^W 2 ^ 4 log q 2 /^ 2 + ■■■. (4.27) 

The key now is to understand precisely which field theory correlator corresponds to taking the 
second derivative of this expression with respect to <p(e)- The field theory action is — -A — F 2 ; 

effectively, we are adding a source by taking the coefficient to be instead — — (1 + Se iq ' x ). 
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Comparing this to Eq. 14.261 we see that 5 = — y ^(e)- Now, we have the two-point correlator 
for G 2 = gy 2 M F 2 : 

J d^(\G 2 (x)\G 2 (0))e^ = ~ {J ^ 1 <i A ^ + • • • , (4.28) 

which should correspond to taking a second derivative with respect to 5 of our above result, 
and so we find: 

.R 3 64N 2 -1 (iV 2 -l) 



k 2 3 64tt 2 3tt 2 
In particular, for N c = 3, this means that 



>5D 



3vr 2 ,2 4 ' 



(4.29) 



(4.30) 



In order to find the actual condensate, we have to differentiate the action with respect 
to the value of the source on the boundary. In our case the source is just the QCD coupling 
itself gy 2 M- Using the expression for the coupling in a pure YM theory 

2n 

- (4 - 31) 

and identifying Aq C d = — we find that the derivative with respect to gy 2 M (viewing g YM as 

a function of z ) is the same as yj^r^o^- 
Putting all this together, we find that: 

Using our estimate of z from the glueball mass in Eq. 14.141 we obtain: 

/^TrF 2 ) = il.19 x 680 MeV) 4 w 0.043 GeV 4 . (4.33) 



\47r 2 



7T 2 



For comparison, an SU(3) lattice calculation found (^G 2 ) w 0.10 GeV 4 [26]. Our result 
is of the same order but slightly smaller. Most phenomenological estimates are smaller, 
beginning with the SVZ result of (^C 2 ) ~ 0.012 GeV , but for pure Yang-Mills the value 
is expected to increase [1]. 



4.3 Relation to Analytic Perturbation Theory 

The QCD perturbation series is an asymptotic expansion of some unknown function, and 
the divergence at Aqcd signals only a breakdown of perturbation theory, not a meaningful 
infinity. In particular, it has been proposed that the pole of the logarithm be cancelled by 
additional terms to produce an "analytic perturbation theory." See, for instance, the work 
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of Shirkov and Solovtsov [27] and related literature (of which there is too much to give an 
exhaustive account here). It is interesting that our holographic equations produce a result 
along these lines, when interpreted in a particular way. 

To see this, note that our identification of the coordinate z with the inverse of a renor- 
malization group scale \i is only clearly defined in the far UV (near z — 0). In fact, when 
we take as a metric ansatz ds 2 = exp(— 2A(y))dx 2 + dy 2 , it is more reasonable to interpret 
A(y) as — log/ii?, so that the 4D part of the metric goes like fi 2 dx 2 . Our expression for the 
QCD coupling in terms of the dilaton exp ^/2/30(y) is given by — ^— , which blows up at a 
finite coordinate in exactly the way that the one-loop QCD beta function tells us a s should 
blow up at /i — Aqcd- On the other hand, using the modified identification of the energy 
scale explained above, 

My) ^ - l °gv R , 

that is by identifying the warp factor (instead of y) as the logarithm of the energy scale, we 
find that y — > yo corresponds to [i — > 0. Thus we can view exp yj2/Z<j){y) as providing a 
formula for a s (fi) that is smoothly defined at all //, which blows up as a power law in the 
deep infrared fj, — > (instead of // — > Aqcd) and reduces to the perturbative result at large 
fi. 

In particular, one can solve for a s (/x) according to this prescription. The relevant expres- 
sions 

*:\») = (y -v(v))/R (4-34) 



lliV, 



c 



can be inverted to find 



-togMi2 = ^ + ^log— (4.35) 
R 4 y - y(fjt) 

Aqcd = z, 1 = ^e^ R (4.36) 



1 llNc -W(^/A QCD ), (4.37) 



where W(y) is the Lambert W-function [28], that is, the principal value of the solution to 
y = xexp(x). In fact, the Lambert W-function appears similarly in the analytic perturbation 
theory approach [29], although the form of a s (fj,) is slightly different there. Nonetheless, our 
results are suggestive of a role for the backreaction on the metric as enforcing good analytic 
properties, which deserves further attention. 



5 Effects of the Tr(F 3 ) condensate 

In pure Yang-Mills, there is an operator of dimension 6, Oq = f abc F®F^ Fp . This operator 
will get a condensate which modifies the OPE at higher orders. We wish to investigate the 
size of the corrections on the glueball masses in our approach, to understand how stable the 
numerics are. Thus we add a new field \i an d we wish to modify the superpotential. We 
should still have in the potential terms A and C exp(60), as before. However, now we should 
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also have a mass term for X , with m 2 x = 6(6 — 4) = 12 by the AdS/CFT correspondence. On 
the other hand, our theory is not conformal, and Oq has an anomalous dimension proportional 
to a s , suggesting we should also have terms in the potential that couple <fi and X . 

For now, we will make no attempt to constrain all the higher-order terms in the 5D action 
coupling and X . Instead, we seek a superpotential with the properties discussed above, as a 
first approximation. Luckily, there is a superpotential which allows the profiles of the scalars 
and the warp factor to be found analytically. Our motivation for this particular choice is 
that it allows an analytic solution and has the desired properties. On the other hand, it 
resembles certain superpotentials that arise in gauged supergravity [30]: 



cxp 




+ cosh(_Bx). 



(5.1) 



(We will see shortly that B — 1.) The corresponding potential is then: 



v(4>,x) 



dW 



dW 
dx 



- 12W 



-12 - Qe^ + (18B 4 - 12B 2 ) X 2 - 3B 2 e^ X 2 + 0( X 4 ). (5.2) 



We find that <f>(y) is as before, whereas x(y) * s given by x'(y) = 67^7 = QB smh(Bx)- But 
this is nearly identical to the equation we solved to find 4>(y) in the case without log running. 
In particular, this means that 



X(y) = log tanh 



3(2 /1 - y) 

R 



(5.3) 



where we have chosen B — 1 to ensure that at small z, X (z) ~ z % . In fact we can check this, 
as in the potential above we expect (18B A — 12B 2 ) X 2 = 6 X 2 , confirming that we want B — 1. 
Finally, we evaluate the warp factor, using A'(y) = W((f)(y), X (y)): 



Mv) = l°g (cosh 3<kV1 r ^ sinh 3<kVi - — ^ + J log 



R 



yo-y 



+ 



m _ log 2 

R 3 : 



(5.4) 



where the first term replaces the y of our solution without the inclusion of Oq, but deviates 
from it in the infrared. (The constant terms correct for a constant difference between y and 
the first term, in the far UV.) The solution in the z coordinates is given by 



ds 2 



1 z ° 
log — 

z 



,12 \ 3 



1 ) dx fJ/ dx u rj IML , — dz 



<t>(z) = - W^loglog^ 



X (z) = log 



2 6 



1 + 



(5.5) 
(5.6) 
(5.7) 
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At this point the reader should be notice certain issues in our calculation. First, while 
it is true that the potential V((f>,x) couples and x an d includes a term suggestive of an 
anomalous dimension, the solutions for (f)(y) and x(y) themselves are completely decoupled! 
This, however, is not really a concern: the backreaction on the metric feels all of the terms 
in the potential. In other words, it is only through A(y) that the anomalous dimension is 
manifest. If, as suggested earlier, we interpret A(y) as — log//, then </>(//) and %(//) will feel 
the effect of the anomalous dimension. 

Another issue is that we have sacrificed precise agreement with perturbation theory for 
the sake of having a simple, solvable example. We have arranged to get the logarithmic 
running of a s , the proper scaling dimension of Oq, and an anomalous dimension term for 
Oq which is proportional to a s . On the other hand, we have not been careful to match the 
coefficient in this anomalous dimension. Details of the OPE and anomalous dimension for 
this operator can be found in Ref. [31]; eventually one would want a model that matches 
them. Of course, in our discussion of cc s (//) earlier, we also had a second (3 function coefficient 
that did not match. This suggests that our analytically solvable superpotentials, while useful 
for a preliminary study, probably need to be replaced by a more detailed numerical study 
based on a more careful matching of the holographic renormalization group. Nonetheless, 
the disagreement appears only at higher orders of perturbation theory, and we can already 
use our preliminary superpotential W(<f>, x) to get some sense of the stability of spectra 
calculated in holographic models. 

5.1 Gubser's Criterion: Constraining z\/zq 

In our solution, <j>(z) blows up at z = z while x( z ) blows up at z — z\. The space will shut 
off at whichever of these is encountered first. Intuitively it is clear that if the dimension 6 
condensate is to make a relatively small correction to the results we have already obtained 
using only the dimension 4 condensate, we should have z 1 > z , so that x( z ) remains finite 
over the interval where the solution is defined. 

In fact there is a conjecture that will enforce this condition. Namely, Gubser in Ref. [32] 
has proposed that curvature singularities of the type arising in the geometries we are con- 
sidering are allowed only if the scalar potential is bounded above when evaluated on the 
solution. By "allowed" , one should understand that in these cases one expects the singular- 
ity to be resolved in the full string theory; geometries violating the criterion are somehow 
pathological. The conjecture is based on some nontrivial consistency checks involving con- 
siderations of finite temperature and examples from the Coulomb branch in AdS/CFT, but 
it is not proven. In any case we will assume for now that it holds for any geometry that 
can be properly thought of as dual to a field theory. It is clear that our original solution 

involving only satisfies the criterion: in that case we had V{<f>) = -12 - Qe^^ < 0. 
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The case with x is more subtle. We have 



V{<f)(z),x{z)) = -15 -6e^ (z) coshx(» + cosh(2x(*)) 




Clearly as Z\ — > oo we recover the previous solution and Gubser's criterion is satisfied. On 
the other hand, as soon as z\ < z , V begins to attain large positive values as z — > Z\. The 
reason is simple: the function 1 — 4x 12 + x has a zero at x ~ 0.9, so the second term above 
can attain positive values when z ~ Z\, and the denominator will attain arbitrarily small 
values provided the singularity at Z\ is reached (i.e. Z\ < z ). In fact large positive values of 
V are attained if and only if z\ < z . 

In short, Gubser's criterion limits us to precisely those solutions which can be viewed in 
some sense as a perturbation of our existing solution. 



5.2 Condensates 



To calculate the condensate we need to again evaluate the classical action for the solution, 
which in our case is given by 



S 



12X| 



•J2 



(5.9) 



+ 12 + 6e V 3<p cosh x + 6 sinlr x ~ K 
Again dropping the UV divergent terms we find either using ()4.19j) or by direct integration 

1 {z\ 2 ^ " 



S 



(5.10) 



2k 2 244 

Following the steps for calculating the condensate for the single field case we find that the 
modified condensate is given by 

(iV 6 2 - 1) 24tt 2 d 



-TrF 2 
4 



3tt 2 llN r Z °'d Zl 



,12 



zl 2 )i 



;i.i9z, 
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(5.11) 



One can see that for z\ > zq (as expected from the criterion of Sec. 15.1)1 this condensate 
is very insensitive to the actual value of Z\. In order to actually fix z\ one would have to 
calculate the second condensate (TrF 3 ) and compare it to the lattice results. However, there 
are no reliable lattice estimates for this condensate available. 



18 



It is also plausible that if we account more properly for the anomalous dimension of 
TrF 3 and for matching of perturbative corrections to the OPE, we will be able to select a 
solution without this ambiguity. However, constructing such a solution appears to require a 
numerical study, which we will leave for future work. 



5.3 Glueball spectra 

Now that we have the background deformed by condensates of Tr F 2 and Tr F 3 , we can again 
compute the glueball spectrum. What we would like to check is how sensitive the results are 
to the value of z\. To simplify the numerical problem we are assuming that the low- lying 
glueball modes are still predominantly contained in the <fi and A fields, and that the leading 
effect of turning on the x field is to modify the gravitational background A(y). Using this 
approximation we find the following equation satisfied by the glueball wave functions (again 
in z coordinates and units of R): 



^4 1 + m 2 z 2 96^ 12 ^ 



F(z) = 0. (5.12) 



One can see that the equation reduces to ()4.13|) in the limit when z\ ^> z > z. By again 
numerically solving this equation for various values of z\/zq > 1 we find that the glueball 
eigenvalues are very insensitive to the actual value of Z\ as long as z\ is not extremely close 
to Zq. For example, the lightest eigenvalue at 2.52/z increases by less than a percent while 
lowering z±/zo from oo to 1.5. For z±/zo = 11 the lightest mass grows by 3 percent, while 
for the extreme value of Zi/zq = 1.01 the growth is still just 9 percent. The glueball mass 
ratios are even less sensitive to Z\\ the ratio of the first excited state to the lightest modes 
decreases by about 3 percent when changing z\/zq from oo to 1.01. Thus we conclude that 
the predictions for glueball spectra presented in the previous section are quite robust against 
corrections from higher condensates. 



6 Linearly confining backgrounds? 

One of the most problematic aspects of holographic QCD is the deep IR physics: one expects 
Regge behavior from states of high angular momentum, and a linear confining potential. The 
solutions presented in this paper show many qualitative and quantitative similarities with 
ordinary QCD. However, they do not properly describe the highly excited glueball states. 
Since there is a singularity at a finite distance, the characteristic mass relation for highly 
excited glueballs will be that of ordinary KK theories (in this respect the theory is similar 
to the models with an IR brane put in by hand) rr? n ~ n 2 /z 2 , instead of the expected Regge- 
type behavior m 2 ~ n/z 2 [8]. Experimental and lattice data suggest that linear confinement 
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effects persist further into the UV than one might expect, and already the light resonances 
observed in QCD seem to fall on Regge trajectories. Regge physics arises naturally from 
strings; in our approach, the more massive excitations of the 5D string correspond to high- 
dimension operators on the field theory side. To describe linear confinement and Regge 
physics accurately, then, it is conceivable that one must take into account the effects of a 
large number of operators. Thus we are led to seek alternative, but still well-motivated, 
approaches to the deep IR physics. One approach is simply to demand that the 5D fields 
have IR profiles that provide the desired behavior, as in Ref. [9]. However, one would like 
to have a dynamical model of this effect. Here we first check that the backgrounds used in 
sections 4-5 do not reproduce such IR profiles. However in [9] it was suggested suggested 
that tachyon condensation might provide the appropriate dynamics. We provide a simple 
modification of our model possibly substantiating this idea, and speculate on its relation to 
known gauge theory effects: namely, UV renormalons and other 1/Q 2 corrections as discussed 
by Zakharov and others [35]. 

6.1 No linear confinement in the dilaton-graviton system 

We first try to find other solutions that would not have a singularity at a finite distance, even 
for the action (|4.1jl since until now we have found only a particular solution for (|4.1|) using 
a superpotential. However in general there should be an infinite family of superpotentials 
giving the same potential. For the general case we cannot find the other solutions analytically 
for all values of the coordinates. We can, however, attempt to solve the equations of motion 
analytically close to the UV boundary. A similar approach has been taken in ref. [33] for 
a type string theory containing a tachyon, considered to be dual to non-supersymmetric 
SU(N) Yang-Mills with six adjoint scalars. In our case we will not consider a tachyon for now, 
as we consider 5D fields to be in one-to-one correspondence with gauge invariant operators 
in the 4D dual. There are some orientifold theories in type that have no bulk tachyon, so 
our approach is not a priori meaningless. 

In order to avoid the change of variables needed to achieve an asymptotically AdS metric 
(plus power corrections), we solve the equations beginning from the explicitly conformally 
flat parametrization of the metric, 



Then for the Einstein equations and the equation of motion for we find (using ' to denote 
derivatives with respect to z): 



ds 2 = e~ 2A{z) (Tj^dzW - dz 2 ) . 






(6.3) 
(6.4) 



3A'(z) 2 - 3A"(z) + (l/2)0'(z) 2 + e- 2Aiz) V{<P) 



0. 



We again use the same potential V(4>) = — 6eV^^( z ) — 6 (but no longer assume the simple 
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expression for the superpotential), and solve the equations in the UV to obtain 



z 



2 log I 8(logi) Z (logi)' 



The omitted terms are those which are smaller as z — > 0, i.e. power corrections in z or higher 
powers in (— logz) -1 . Here 7 is a parameter that is undetermined by the UV equations of 
motion, just as in the solution of ref. [33]. Further, we find that if we add small perturbations 
to exp( v /2730) and to A(z), an 0(z 4 ) power correction to the former is allowed while only 
constant corrections to the latter are allowed. It is reassuring that we find a power correction 
suitable for the gluon condensate. Ref. [33] found also a correction of order z 2 and interpreted 
it as a UV renormalon 1 (whereas the power correction we have considered may be thought 
of as an IR renormalon). The tachyon was crucial to the appearance of the UV renormalon 
term, which might have interesting implications. 

In order to find a linearly confining solution, we would need to connect these UV solutions 
to solutions in the z — > 00 IR region, which would also give linear confinement. Suppose we 
then search for a solution that is valid at large z. In the case that there is no cosmological 
constant, there is a linear dilaton solution; this is clear since our equations are those for a 
noncritical string. The linear dilaton persists in the presence of a cosmological constant as a 
solution at large z; the equations are satisfied up to terms exponentially small at large z by 
a linear dilaton on flat 5D space. Such a linear dilaton background does not give a confining 
solution. Furthermore, one can check that no other power-law growth in z is allowed. In 
particular, the z 2 behavior as in Ref. [9] is not a solution to our equations of motion. Thus 
we conclude that none of the solutions of the action ()4.1|) will result in a linearly confining 
background. 



6.2 Linear confinement from the tachyon-dilaton-graviton system? 

Although we do not understand how to systematically approach the study of UV renormalon 
corrections, it is at least superficially plausible that they are associated with a closed string 
tachyon. This is a natural scenario to consider, since the UV renormalon is associated with 
Q~ 2 corrections. Ref. [33] found such an association in a concrete string theory model. The 
study of such corrections has been discussed in great detail in the QCD literature, which 
links the idea of the UV renormalon to the quadratic corrections associated with the QCD 
string tension and with a hypothetical nonperturbative tachyonic gluon mass, as well as 
with monopoles and vortices; we can only refer to a small sample of that literature [35]. 
Holographically, effects associated with dimension 2 corrections would appear to be associ- 
ated with tachyonic physics. Interestingly, it has been proposed that the UV renormalon 
is associated with a nonvanishing value of (A 2 ) min where the minimization is over gauge 

1 For a review of the physics of the UV renormalon, see Ref. [34]. 
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choices [36]. Such a dimension 2 condensate could plausibly be associated holographically to 
a closed string tachyon (one that saturates the Breitenlo liner- Freedman bound [37]), though 
it does contradict the picture of holographic fields as being associated to gauge-invariant, 
local operators in the field theory. The minimum value of the A 2 condensate can be ex- 
pressed in terms of nonlocal, gauge-invariant operators, beginning with F IXU (D 2 )~ 1 F^ v [38]. 
The nonlocality of this operator is perhaps suggestive of the long-distance, stringy effects of 
the flux tube needed to describe excited hadrons. 

All of these hints of the importance of Q~ 2 corrections suggest that we take the idea of a 
closed string tachyon seriously, despite the lack of a gauge-invariant local operator for it to 
match to. Perhaps this indicates that certain 5D corrections can somehow be "resummed" 
into the effects of a tachyonic field. The linear confining potential of QCD has been suggested 
to relate to an expcz 2 background in 5D [9]. We would like to have a dynamical solution 
incorporating both asymptotic freedom in the UV and linear confinement in the IR. We have 
observed that the cosmological constant does not destroy the existence of a linear dilaton 
solution at large z (up to small corrections). This suggests that if we have a theory with a 
quadratic tachyon profile at zero cosmological constant, such a solution may persist in the 
presence of a cosmological constant. 

The action of the bosonic noncritical string theory including the leading a 1 corrections, 
using a sigma model approach, is (after transforming to Einstein frame), according to 



4 $ 4 



2TZ 



a 



|(34) 2 



a 



--(<9$) 2 - -d$dt+ {dtf 



(6.7) 



To this action we are adding a cosmological constant (which we assume could come from a 
5- form flux) and adjust the parameters such that with the t — > substitution we recover the 
action considered in the previous section. The resulting action is: 
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2^2 



d 5 Xy/g 



4 $ 12A-1 12 
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(6.8) 



The resulting equations of motion for the metric ansatz ds 2 = e 2A<yZ \dx 2 — dz 2 ) are: 



3 R 2 \ 



(A - 1) + f(3A>& + $") + e -« (^* +2A |^(l + t~ 7^)^ + 



(6.9) 



12 



- - t) - 12A' 2 - 8A" - -$' 2 - t' 2 - 4A'& - + 2t" + 6A't' = 



(6.10) 
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-6A' 2 (2 - e-*) + ±<M(1 + l e f(A - 1 + e-*(l - ^ + t)) + $' 2 (jj - ^) + 
1 -e- t (-U't' + t' 2 ) = 0. (6.11) 

Here A = (26 — D)/3 = 21/3. One can show that for large z the leading order solution to 
these equations is given by 

$^$ , A^A , t ^-L 2A ° + l*°z 2 + l(\-2). (6.12) 

A 2 

In order to serve our purposes there must be a solution interpolating between logarithmic 
running on AdS in the UV and the above confining background (with flat space and constant 
dilaton) in the IR. It would be interesting to numerically study these equations, as well as 
similar systems in the type string. The IR solution above has the property that the e~* 
terms are growing at large z, whereas the proposed background in Ref. [9] has an exponential 
shutoff of the metric, which is not a solution of the above equations. However, the sign does 
not appear to affect the existence of Regge physics. Regardless, we intend these remarks not 
as a definitive statement, but as a provocative hint that further studies of tachyon dynamics 
could be useful for understanding QCD. 



7 Conclusions and Outlook 

In this paper we have clarified a number of aspects of AdS/QCD, and sketched a concrete 
program for computing in the holographic model and estimating associated errors. We have 
also explained that deep IR physics, associated with high radial excitations or large angular 
momentum, is troublesome in this framework. The underlying reason is that the OPE does 
not reproduce the quadratic corrections associated with this physics. We have suggested that 
a closed string tachyon can reproduce much of the underlying dynamics, but at this point 
that is a toy model and it is not clear how to systematically apply the idea to computations. 

Our results suggest a number of directions for new work. One obvious task is to extend 
these results to theories with flavor. This should be straightforward, although there are 
potential numerical difficulties. It would be particularly interesting to see if one can obtain 
results for mixing of glueballs with qq mesons without large associated uncertainties. 

Another direction is to make more explicit the connection between the renormalization 
group and the holographic dual. The 5D action has an interpretation as the 4D generating 
functional W[J], whereas some numerical attempts involving truncations of exact RGEs 
have focused on computing the Legendre transform of this quantity, There are subtle 

issues of nonperturbative gauge-invariant regulators that must be considered in such studies, 
but it is conceivable that such existing work could be reinterpreted as a computation of 
background profiles for various fields, about which we could then compute the spectrum of 
excitations and couplings with the usual 5D techniques. We hope to clarify this relationship 
in a future paper. It is also interesting in this context to think about how the holographic 
renormalization group might relate to the "analytic perturbation theory" framework. 
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Finally, the intricate story relating 1/q 2 corrections, UV renormalons, vortices, linear 
confinement, and the closed string tachyon is very interesting and still rather poorly under- 
stood. A better understanding of these quantities and their relationships could elucidate the 
confining dynamics of QCD, and also shed light on closed string tachyon condensation in 
general (with possible applications to cosmology). 
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